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ABSTRACT When challenged with osmotic shock, Saccharomyces cerevisiae induces hundreds of genes, KEYWORDS 
despite a concurrent reduction in overall transcriptional capacity. The stress-responsive MAP kinase Hog1 MAP kinase 
activates expression of specific genes through interactions with chromatin remodeling enzymes, transcription transcription 
factors, and RNA polymerase II. However, it is not clear whether Hog1 is involved more globally in modulating stress response 
the cell's transcriptional program during stress, in addition to activating specific genes. Here we show that osmotic shock 
large-scale redistribution of RNA Pol II from housekeeping to stress genes requires Hog1. We demonstrate Hog1 
that decreased RNA Pol II occupancy is the default outcome for highly expressed genes upon stress and that 
Hog1 is partially required for this effect. We find that Hog1 and RNA Pol II colocalize to open reading frames 
that bypass global transcriptional repression. These activation targets are specified by promoter binding of two 
osmotic stress-responsive transcription factors. The combination of reduced global transcription with a gene- 
specific override mechanism allows cells to rapidly switch their transcriptional program in response to stress. 



Saccharomyces cerevisiae responds to many environmental perturba- 
tions with large-scale changes to its transcriptional program (Gasch 
et al. 2000). For example, when yeast cells are subjected to osmotic 
shock, the resulting transcriptional response rapidly alters mRNA 
levels of over one third of the genome (O'Rourke and Herskowitz 
2004). During the initial phase of the osmotic shock response, the 
global transcription rate drops by 50% (Romero-Santacreu et al. 2009) 
due to a general defect in transcription initiation in the first few min- 
utes of stress (Proft and Struhl 2004). Despite this global reduction, 
transcription of hundreds of genes increases, and the most strongly 
activated genes are induced over 100-fold within 20 min (Capaldi 
et al. 2008; O'Rourke and Herskowitz 2004), suggesting that cells 
may have evolved a strategy to prioritize transcription of these genes 
in the face of harsh conditions. 
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Transcriptional activation of osmotic stress response genes is co- 
ordinated by the MAP kinase Hogl (Capaldi et al. 2008; O'Rourke 
and Herskowitz 2004). Upon stress, Hogl is localized to the nucleus 
(Ferrigno et al. 1998) where it regulates the action of several transcription 
factors (Alepuz et al. 2001; Proft et al. 2001; Rep et al. 2000) and 
chromatin-modifying enzymes (De Nadal et al. 2004; Mas et al. 2009). 
In addition, previous work strongly suggests that Hogl interacts directly 
with RNA Pol II in stressed cells. This interaction occurs both in vitro 
(Alepuz et al. 2003), as purified Hogl and RNA Pol II interact in the 
absence of DNA, and in vivo (Alepuz et al. 2001; 2003; Proft et al. 2006), 
even in cells lacking Hogl kinase activity, or when stress-responsive 
transcription is blocked by deletion of stress responsive transcription 
factors (Alepuz et al. 2003). In stressed cells, Hogl interacts with com- 
ponents specific to the RNA Pol II initiation complex (Alepuz et al. 

2003) and is sufficient to recruit RNA Pol II when artificially tethered 
to a promoter (Alepuz et al. 2003). In addition to its role in transcription 
initiation, Hogl interacts with the elongating RNA Pol II holoenzyme 
(Proft et al. 2006) and colocalizes with RNA Pol II in the open reading 
frames (ORF) of some genes (Pokholok et al. 2006; Proft et al. 2006). A 
similar behavior has been observed for another stress-induced yeast 
MAP kinase, Mpkl, which acts as a gene- specific elongation factor 
during heat shock by blocking attenuation (Kim and Levin 2011). 
Hogl regulates chromatin state at specific promoters (De Nadal et al. 

2004) , promotes transcription initiation at specific genes (Alepuz et al. 
2003), and may, like Mpkl (Kim and Levin 2011), act as a gene-specific 
elongation factor, which has been suggested previously (Proft et al. 2006). 
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However, it is not clear whether Hogl plays a role in coordinating the 
global changes in transcription that accompany osmotic stress (Miller 
et al 2011; Proft and Struhl 2004; Romero-Santacreu et al 2009), in 
addition to activating transcription of specific genes. 

In this study, we investigate the role of Hogl in global transcrip- 
tional allocation of resources in the early stages of stress. In agreement 
with previous work (Miller et al 2011; Romero-Santacreu et al 2009), 
we observe a global reallocation of RNA Pol II upon osmotic shock, in 
which housekeeping genes lose polymerase occupancy as stress genes 
are activated. We demonstrate that this reallocation of RNA Pol II 
depends upon Hogl. We find that depletion of RNA Pol II is the 
default outcome for highly expressed genes in response to stress and 
that this depletion is less pronounced in the absence of Hogl. We find 
that RNA Pol II and Hogl colocalize to the ORFs of the most highly 
induced genes during stress. Hogl colocalizes with RNA Pol II specif- 
ically at a set of stress-responsive ORFs that are marked by stress- 
induced promoter binding of Hogl and two of its cognate transcription 
factors, Skol and Hotl. We propose a model in which RNA Pol II is 
preferentially recruited to promoters that contain Skol and Hotl bind- 
ing sites, allowing recruitment or assembly of a Hogl -RNA Pol II 
complex and prioritizing transcription of these genes during stress. 

MATERIALS AND METHODS 
Strains 

All strains used in this study are in the W303 strain background (trpl 
leu2 ura3 his3 canl GAL+ psi+), and are listed in supporting infor- 
mation, Table SI. Epitope tags and gene deletions were introduced by 
transformation with PCR products including auxotrophic or antibiotic 
markers flanked by the 40 base pairs (bp) of sequence found directly 
upstream and downstream from the gene, followed by selection on the 
appropriate medium (Longtine et al 1998; Rothstein 1991). Strains 
with multiple gene manipulations were constructed by mating the 
single deletion strains and dissecting the resulting tetrads. 

A Quickchange kit (Stratagene) was used introduce amino acid 
substitutions to generate nonphosphorylatable Skol variants. First, 
3HA-Skol was cloned out of genomic DNA into YCp50, where muta- 
tions were introduced. Mutant Skol was then reintroduced at the 
endogenous locus. 

Stress-promoted LACZ constructs were derived from pCM173 
plasmid (Gari et al 1997), obtained from EUROSCARF. The region 
between the EcoRI site and the 5' end of LACZ was replaced with 
a promoter region of interest. 

Chromatin immunoprecipitation 

ChlP-seq and ChlP-qPCR experiments were conducted as described 
previously (Capaldi et al 2008; Johnson et al 2007; Robertson et al 
2007), and methods are described in detail in the supporting information. 
Briefly, cells were grown at 30° in YEPD (or SD for experiments re- 
quiring selection for a plasmid), with shaking, to OD 600 of 0.6. For stress 
treatment, YEPD supplemented with KC1 was added to a final concen- 
tration of 0.4 M; for mock-treated cells, the same volume of YEPD was 
added to cultures. Samples were harvested after 5 min of stress or mock 
treatment. ChIP against HA-tagged strains was performed with the 
monoclonal anti-HA antibody 12CA5, and for RNA Pol II ChIP, 1Y26 
anti-Rpb3 (Neoclone) antibody was used. For ChlP-seq, ~10 ng IP 
material was used to generate each library, following the Illumina pro- 
tocol for their paired end DNA sample prep kit (vl). For ChlP-qPCR, 
samples were analyzed on an MX3000p qPCR machine (Stratagene) 
using primers that amplify a ~100 bp region surrounding the center 
of observed binding peaks (primer sequences are listed in Table S2). 



More detailed methods and analysis of ChIP data are presented in the 
supporting information. ChlP-seq results have been confirmed by qPCR 
at selected genes (PDC1, JLV5, TDH3, ADH1, PMA1, STL1, RTC3, 
HSP12). Raw data for all ChlP-seq experiments may be obtained from 
Gene Expression Omnibus (GEO) under accession number GSE38208. 

Inducible LACZ 

YCP50 plasmid carrying P TET07 -IACZ, cloned from pCM173 plas- 
mid (Gari et al 1997), obtained from EUROSCARF, was produced 
and transformed into a strain expressing rtTA from the MY02 pro- 
moter, integrated at the URA3 locus (Kim and O'Shea 2008). Doxy- 
cycline concentrations ranging from 0.1 to 2 jig/mL -1 were tested to 
determine an appropriate working concentration that would drive 
RNA Pol II occupancy of the LACZ construct to a level similar to 
that of the most highly expressed endogenous genes. Based on ap- 
proximately 5 hr of growth in rich medium at 30° after addition of 
doxycycline, 1.25 u>g/ mL" 1 doxycycline was determined to be an 
appropriate concentration of inducer. For the experiment, cells were 
grown to saturation overnight, and then diluted to OD 600 0.1. At this 
point, doxycycline was added to induce LACZ expression. 

Hogl -occupied ORFs 

We determined Hogl enrichment at each ORF by summing Hogl 
ChIP reads that align anywhere in the ORF, subtracting reads 
obtained in a mock IP, and dividing by reads from sequenced input 
material. We used three criteria to designate Hogl ORF occupied 
genes: Hogl enrichment; ratio of Hogl ChIP signal in the presence 
and absence of stress; and length-normalized Hogl ChIP signal within 
the ORF. Cutoff for each criterion was three standard deviations above 
the genome-wide median. To estimate the false-positive rate, we ap- 
plied these same criteria to Hogl ChlP-seq data in the absence of 
stress (when Hogl is largely excluded from the nucleus); we identify 
only one gene that passes all three criteria. We excluded the 5% of 
genes with the lowest number of input reads to avoid variation in 
signal due to low alignability. Where two ORFs overlap, we use the 
nonoverlapping portions to calculate separate enrichment values. 
Nineteen genes on our list overlap with those identified previously 
by a lower resolution ChlP-chip assay (Pokholok et al 2006). 

ChlP-seq peak identification 

Initial identification of peaks was performed using an implementation 
of the PeakSeq method (Rozowsky et al 2009) written by Xu Zhou 
(Zhou and O'Shea 2011). This list of peaks was screened by additional 
criteria: for one hundred consecutive bases, peaks must be two input 
standard deviations about the input value at that location and above 
the overall median input value. This screened list of peaks was then 
sorted by enrichment over input, and the one hundred highest enrich- 
ment peaks were selected for analysis. For this purpose, enrichment 
was denned as the ratio of IP reads to input reads, with reads summed 
over the 50 base pairs at the center of each binding peak. Complete lists 
of all peaks found using PeakSeq are available from Gene Expression 
Omnibus (GEO) under accession number GSE38208. 

Identification of Skol and Hotl candidate motifs 

We conducted a bioinformatics search of the genomic regions enriched 
in each ChlP-seq experiment using the MEME (Bailey and Elkan 1994) 
and MAST (Bailey and Gribskov 1998) tools of the web-hosted MEME 
software suite (Bailey et al 2009). We selected 50 bp regions surround- 
ing the maximum height position of each peak for analysis. During 
motif discovery, we specified the parameter that each sequence should 
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Figure 1 Reallocation of RNA Pol II upon osmotic shock. 
(A) Dynamics of depletion and recruitment of RNA Pol II in 
stress were measured by a time course of ChlP-QPCR 
experiments. Cells in log growth were subjected to 
osmotic shock by addition of 0.4 M KCI to the growth 
medium, and samples were collected at the indicated time 
points. Housekeeping genes with high RNA Pol II occu- 
pancy pre-stress are quickly depleted of Pol II upon stress 
and begin recovering within 10 min (average behavior of 
PDC1, TDH3, ILV5). Stress-induced genes recruit polymer- 
ase upon stress treatment and then rapidly return to low 
RNA Pol II occupancy (average behavior of RTC3 and 
STL 7). Data shown are the averages of biological replicate 
experiments. (B) RNA Pol II occupancy was measured by 
ChlP-sequencing of Rpb3, the third largest subunit of RNA 
Pol II, in control (mock-treated with YEPD) cells and in cells 
subjected to 5 min of osmotic shock induced by addition 
of 0.4 M KCI to the growth medium for 5 min. Histograms 
show the distribution of RNA Pol II occupancy across all 
yeast ORFs in each condition. RNA Pol II occupancy is 
measured at each ORF in reads per kilobase per million 
(RPKM), after subtracting the input value. Regions of over- 
lap between two or more ORFs were excluded from anal- 
ysis. Inset plots show the 95 th percentile for RNA Pol II 
occupancy in each condition; the 95 th percentile is indi- 
cated by a dashed line in the main plot. (C) Scatter plot 
showing RNA Pol II occupancy during rapid growth in 
YEPD (x-axis) vs. osmotic shock in 0.4 M KCI for 5 min 
(y-axis). The number of reads that align to each ORF is 
normalized to RPKM. Each point represents one gene, 
and points are color-coded to show ORFs that are occu- 
pied by Hog1 in response to stress (red points). (D) Distinct 
sets of genes are top RNA Pol II targets in the presence 
and absence of stress. RNA Pol II occupancy in rapid 
growth (gray plots) and after osmotic shock for 5 min in 
0.4M KCI (black plots). Values plotted are the number of 
reads that align at each point along the chromosome, and 
each ChlP-seq dataset is scaled to one million reads. 
Arrows indicate the position and direction of each open 
reading frame. The five highest occupancy ORFs during 
osmotic shock are shown. (E) The five highest occupancy 
ORFs during rapid growth in YEPD are plotted, as in (D). 



contain one occurrence of the putative binding motif. We scanned 
peaks that are bound by Hotl in stress, bound by Skol with higher 
enrichment pre-stress, and bound by Skol with higher enrichment in 
stress. For discovery of the stress-induced Skol binding motif, peaks 
bound by Skol in the absence of stress were used as a set of counter- 
examples in a discriminatory search. Matches to the discovered Skol 
and Hotl motifs in the regulatory regions surrounding Hog 1 -occupied 
ORFs are listed in File S2 Binding data and motifs found for Skol and 
Hotl are in general agreement with previous measurements made by 
ChlP-chip under similar conditions (Capaldi et al 2008). 

RESULTS 

Global reallocation of RNA Pol II during osmotic shock 

To investigate allocation of transcriptional resources during the early 
stages of the S. cerevisiae osmotic stress response, we measured RNA 
Pol II location genome-wide by chromatin immunoprecipitation read- 
out by high-throughput sequencing (ChlP-seq) of Rpb3, the third larg- 
est subunit of RNA Pol II. We assume that RNA Pol II ChIP signal is 
proportional to the instantaneous transcription rate, as has been dem- 



onstrated previously (Miller et al 2011). Osmotic shock was induced by 
addition of 0.4 M KCI to the growth medium, and samples were col- 
lected after 5 min of stress treatment. We selected this time point to 
observe the transition from the initial shock phase of the stress response 
to the induction phase — at 5 min post-stress, the global transcription 
rate is reduced, but stress responsive genes are actively transcribed 
(Miller et al 2011; Romero-Santacreu et al 2009; and Figure 1A). Pre- 
vious work (Miller et al 201 1) used dynamic labeling of nascent RNA to 
characterize a redistribution of RNA Pol II upon osmotic shock that 
alters transcriptional activity of about 10% of the genome. Using a high- 
resolution measurement of RNA Pol II occupancy, we observe changes 
in RNA Pol II occupancy that are consistent with these observations 
(Figure SI), and additionally, we note that while these stress-induced 
changes in transcriptional activity affect a relatively small fraction of 
genes (Figure 1, B and C), this redistribution of RNA Pol II constitutes 
a major reallocation of the cell's transcriptional resources (Figure 1C). 
After 5 min in stress, a set of stress-induced genes shows the highest 
RNA Pol II occupancy out of all genes in the genome (Figure ID,), 
while RNA Pol II occupancy decreases substantially at the genes that 
had been highly transcribed (Figure IE). 



:~:£G3'Genes | Genomes | Genetics 



Volume 2 September 2012 I Reallocation of RNA Pol II upon Stress I 1131 



RNA Pol II is depleted from a heterologous gene upon 
osmotic shock 

Why are high expression genes depleted of RNA Pol II upon stress, 
while stress genes are simultaneously induced to a very high level? 
Perhaps depletion of RNA Pol II is the default response to stress, due 
in part to the physical effects of stress (Proft and Struhl 2004). To ask 
if stress-induced depletion of RNA Pol II from highly expressed genes 
requires specific ds-elements, we drove the expression of an exoge- 
nous, plasmid-encoded gene (LACZ) with an exogenous induction 
system (Tet-On) and asked if this gene would lose RNA Pol II occu- 
pancy in response to stress. We then measured RNA Pol II occupancy 
in the presence or absence of osmotic shock (5 min in 0.4 M KC1) at 
LACZ and at two control genes (Figure 2A). As expected, RNA Pol II 
is depleted from the housekeeping gene ADH1 upon stress, and 
recruited to the stress-induced gene RTC3. Induced LACZ behaves 
like the endogenous housekeeping gene control: RNA Pol II is highly 
enriched at LACZ before stress, and is depleted from LACZ upon 
stress. This result suggests that the default response to stress is de- 
pletion of RNA Pol II, and that the coding or regulatory regions of 
highly transcribed stress-induced genes contain a'5-elements that en- 
able them to bypass this effect. 

Hog1 controls redistribution of RNA Pol II upon stress 

In cells lacking Hogl, the decrease in total transcription that occurs 
within minutes of stress exposure is initially less pronounced 
(Romero-Santacreu et al 2009), suggesting a possible role for Hogl 
in global transcriptional regulation. To ask if Hogl is required for 
depletion of RNA Pol II from highly expressed genes upon stress, 
we tested the effect of osmotic shock on expression of our heterolo- 
gous LACZ expression system in a strain lacking Hogl (Figure 2A, red 
bars). In Ahogl cells, a greater proportion of RNA Pol II is retained at 
the housekeeping gene ADH1 upon stress treatment, and as expected 
(Capaldi et al 2008), there is a defect in Pol II recruitment to the Hogl 
regulated stress-responsive gene RTC3. As in wild- type cells, the het- 
erologous LACZ reporter behaves like the endogenous housekeeping 
gene control, showing less severe depletion of RNA Pol II upon stress 
in the absence of Hogl than in wild-type cells. Thus, depletion of 
RNA Pol II from genes upon stress is partially dependent on Hogl. 

To test the influence of Hogl on the global reallocation of RNA 
Pol II that we observe in stress (Figure 1C), we measured RNA Pol II 
occupancy genome-wide in the presence and absence of stress in 
a Ahogl strain. In the absence of Hogl, we no longer observe re- 
distribution of RNA Pol II upon stress (Figure 2B and Figure S2). 
Thus, the decrease in the overall transcription rate affects all genes 
similarly in Ahogl cells, while wild- type cells actively prioritize tran- 
scription of stress-responsive genes over housekeeping genes. 

Under osmotic shock conditions, Hogl is found in a complex with 
RNA Pol II (Alepuz et al 2003; Proft et al 2006) and is present in 
some ORFs during elongation (Pokholok et al 2006; Proft et al 2006). 
To investigate the relationship between high RNA Pol II occupancy 
and the presence of Hogl in ORFs during stress, we determined Hogl 
localization genome- wide by ChlP-seq in the presence and absence of 
osmotic stress. We identified 28 ORFs that are enriched for Hogl 
during osmotic shock (listed in File SI). These findings are in general 
agreement with previous work (Pascual-Ahuir et al 2006; Pokholok 
et al 2006; Figure S3), but by using a higher resolution technique, we 
are able to more clearly differentiate between Hogl present in regu- 
latory regions vs. coding regions. We find that Hogl occupies the 
ORFs that show the highest increase in RNA Pol II occupancy upon 
stress (Figure 1C, red points), while Hogl is not present in ORFs that 
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■ wtKCI in YEPD, RPKM 

■ Ahogl KCI 

Figure 2 Reallocation of RNA Pol II upon osmotic shock is Hog1- 
dependent. (A) RNA Pol II is depleted from a heterologous gene upon 
stress. The TET-ON system was used to drive expression of Pteto7 _ 
LACZ, encoded on a plasmid, to a high level in wild-type cells or in 
a strain lacking Hogl. Induced cells were then subjected to osmotic 
shock (as described in (A). RNA Pol II enrichment over a featureless 
(subtelomeric) region was measured at LACZ and at two control genes: 
ADH1, which shows high RNA Pol II occupancy in the absence of 
stress, and RTC3, which shows stress-induced RNA Pol II occupancy. 
Values plotted are relative RNA Pol II enrichment values for each gene. 
Error bars show the standard deviation of three biological replicates. 
(B) Scatter plot showing RNA Pol II occupancy in a Ahogl strain during 
rapid growth in YEPD (x-axis) vs. osmotic shock in 0.4 M KCI for 5 min 
(y-axis). Each point represents one gene. The number of reads that 
align to each ORF are normalized to reads per kilobases per million 
(RPKM). 

show a decrease in RNA Pol II occupancy in response to stress or in 
the ORFs of HOG pathway-activated genes that show only modest 
induction upon stress. Hogl -occupied ORFs show higher RNA Pol II 
occupancy (median value of 17.1 -fold enrichment over input) than 
HOG pathway-activated genes that are not occupied by Hogl (median 
value of 2.4-fold enrichment over input). 

Promoter regions are sufficient to direct Hogl to ORFs 

Previous work demonstrates that for at least one Hogl -occupied ORF, 
the region downstream of the ORF is sufficient to induce Hogl ORF 
occupancy (Proft et al 2006). However, it is not clear which cis- 
elements confer this ability, or whether specific transcription factors 
mediate Hogl recognition of these c/s- elements. Hogl localizes to 
promoter regions of some target genes (Alepuz et al 2001, 2003; 
Capaldi et al 2008; Proft and Struhl 2002) as well as ORFs, but the 
relationship between these two events is unclear, and Hogl recruit- 
ment to ORFs is not thought to be associated with any particular 
transcription factor (Proft et al 2006). We find that Hogl presence 
in ORFs correlates strongly with Hogl promoter binding. Hogl shows 
enrichment (denned as mock- subtracted Hogl ChIP signal divided by 
input for the region 1000 bp upstream of each ORF) greater than three 
standard deviations above the genome- wide median in the promoters 
of 26 of the 28 occupied ORFs (example genes shown in Figure 3A), 
compared with 5 of the 231 promoters of HOG pathway- regulated 
ORFs that are not occupied by Hogl. We also observe Hogl ChIP 
enrichment greater than three standard deviations above the genome- 
wide median in the 3' regions (50-500 bp downstream) of 18 of the 
28 ORFs that are occupied by Hogl, compared with 4 out of 231 of 
the HOG-pathway regulated ORFs that are not occupied by Hogl. 

We were interested in whether promoter regions contain in- 
formation sufficient to direct Hogl to ORFs. To address this question, 
we placed a stress gene promoter upstream of LACZ on a plasmid and 
tested for the presence of Hogl in the LACZ ORF. We find that the 
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Figure 3 The Hog1 cognate TFs Sko1 and Hot1 target 
stress-induced genes for transcription by a Hog1-RNA 
Pol II complex. (A) Colocalization of Hog1, Sko1, and 
Hot1 at example genes occupied by Hog1 during 
stress. Lighter colored overlay for Pol II and Hog1 
ChlP-seq shows ChIP signal in the Asko1Ahot1 strain. 
Black arrows and dashed gray guidelines mark the po- 
sition of the GPD1, RHR2, STL1 , and SCM4 ORFs, and 
white triangles show predicted transcription factor bind- 
ing sites. Values plotted are the number of reads that 
align at each point along the chromosome, and each 
ChlP-seq dataset is scaled to one million reads. High- 
lighted regions are those shown to be sufficient for 
Hog1 ORF occupancy. (B) RNA Pol II and Hog1 pres- 
ence in the LACZ ORF driven by three different pro- 
moters (defined as 1000 base pairs upstream of each 
gene): the stress gene promoters ?gpdi and Prhr2, and 
Pcycu an inactive promoter that serves as a negative 
control. These promoter- LACZ constructs are carried 
on a plasmid. Hog1 and RNA Pol II ChIP enrichment 
was measured by qPCR during rapid growth in SD me- 
dium lacking tryptophan to maintain the plasmid (bars 
labeled -) and in response to osmotic shock induced by 
addition of 0.4 M KCI to the medium (bars labeled +). 
Enrichment values are normalized by ChIP to POLL For 
Hog1 enrichment, no-stress values were normalized to 
one, and stress values show the relative change in en- 
richment. Error bars show the standard deviation of 
three biological replicates. 



promoters of the stress genes GPD1 and RHR2 are sufficient to recruit 
both RNA Pol II and Hogl to the LACZ ORF in a stress- specific 
manner (Figure 3B). This result, combined with the observation that 
the region downstream of STL1 is necessary and sufficient to recruit 
Hogl to the STL1 ORF (Proft et al 2006), suggests that cis -regulatory 
elements contained within these 3' and 5' regions are capable of 
targeting genes for transcription by a complex that includes Hogl. 

ORFs occupied by Hogl are marked by the presence 
of Hogl, Sko1, and Hot1 in their regulatory regions 

Although Hogl lacks a DNA binding domain, it is able to associate 
with promoters through physical interactions with its cognate tran- 
scription factors Skol (Proft and Struhl 2002) and Hotl (Alepuz et al. 
2001, 2003). Of the 28 ORFs occupied by Hogl, 21 require Skol and/ 
or Hotl for activation in response to osmotic shock (Capaldi et al. 
2008). We hypothesized that these factors may be able to recruit 
Hogl, and thus the presence of these transcription factors at pro- 
moters would correlate with Hogl recruitment to promoters and 
ORFs. To obtain a high-resolution map of Skol and Hotl binding 
within the regulatory regions surrounding Hogl occupied ORFs, we 
measured Skol and Hotl binding by ChlP-seq in the presence and 
absence of stress. At promoters where Hogl is present, Skol and Hotl 
colocalize with Hogl (Figure 3A). Skol is present in the promoters of 
27 of the 28 Hogl-occupied ORFs, and Hotl binds within 50 bp of 
Skol at 14 of these promoters. At each of these promoters, the ob- 
served transcription factor binding peaks overlap with Hogl promoter 
peaks. 

Using a bioinformatics approach, we identified sequence elements 
that are overrepresented at sites of Hotl and Skol recruitment in 



response to stress (Figure 4A). We find that Hogl colocalizes with 
Skol and Hotl at overrepresented sites in the promoters, and in some 
cases, downstream regions of the ORFs that it occupies (Figure 3A, 
white triangles). These predicted Skol and Hotl binding sites are 
present in regulatory regions that are sufficient to induce Hogl 
ORF occupancy (Figure 3A, highlighted regions). Hogl recruitment 
to these regulatory regions depends upon the presence of Skol and 
Hotl; Hogl is not recruited to promoters, ORFs, or downstream 
regions in the absence of these two factors (Figures 3 A and 4B). 

Hogl is required for stress-induced Skol 
and Hotl binding 

Although Hogl phosphorylates Hotl upon stress, this phosphoryla- 
tion is not required for Hotl activity (Alepuz et al. 2003). Instead, 
a physical interaction between Hotl and Hogl is required for pro- 
moter recognition, and it is Hogl that recruits the polymerase (Alepuz 
et al. 2003). Skol binds constitutively to some genomic locations, 
while Skol binding at other sites is induced by stress (Capaldi et al. 
2008; Ni et al. 2009). While Hogl is not required for Skol binding in 
the absence of stress, we hypothesize that Hogl may play a similar role 
in directing stress-induced binding of Hotl and Skol, given that Hogl 
colocalizes with Skol at stress-induced binding sites (Figure 3A), and 
Skol is not recruited to these sites upon stress in a Ahogl strain 
(Figure 4C). Hogl phosphorylates Skol in response to stress (Proft 
and Struhl 2002), but we find that this phosphorylation is not required 
for stress-induced Skol binding (Figure 4C), suggesting that Hogl 
regulates Skol binding as well as Hotl binding through a sustained 
physical interaction rather than via phosphorylation. Though Hogl 
lacks a DNA binding domain, it partners with Skol and Hotl to 
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Figure 4 Analysis of binding 
determinants of Sko1 and Hot1. 
(A) Binding behavior (left panel) 
and discovered motifs (right 
panel) for Hot1 and two classes 
of Sko1 binding sites observed by 
ChlP-seq. Sko1 and Hot1 binding 
is shown in control conditions and 
in osmotic shock (treatment with 
0.4 M KCI for 5 min. Sko1 no- 
stress peaks show decreased en- 
richment upon osmotic shock 
(example peaks from FSH1 and 
PHD1 promoters, top left panel.) 
Stress-induced Sko1 (center left 
panel) and Hot1 (bottom panel) 
binding peaks increase in ChIP 
enrichment upon stress (example 
peaks in STL1 and RTC3 pro- 
moters). A bioinformatics search 
of 43 pre-stress Sko1 binding 
sequences reproduces a previ- 
ously reported Sko1 binding mo- 
tif, 5'-TGACGT-3' (Harbison etal. 
2004). A bioinformatics search of 
17 stress-induced Sko1 binding 
peak sequences reveals a low- 
information variant of the Sko1 

consensus site. A search of 11 high-confidence Hot1 binding peaks allows de novo discovery of a candidate motif for Hot1. (B) Hog1 enrichment during 
stress at ORFs, promoters, and 3' regions decreases in the absence of Sko1 and Hot1. Box and whisker plot shows Hog1 enrichment in the ORFs, 
promoters (1000 bp upstream) and 3' regions (1000 bp downstream) for the set of 28 ORFs occupied by Hog1 during stress in wild-type cells. Enrichment 
is defined as Hog1 ChlP-seq reads minus mock IP, divided by the input signal. Enrichment values in wild-type (wt) cells are shown in dark blue; and 
enrichment values for Asko1Ahot1 cells are shown in light blue. The edges of each box show the 25 th and 75 th percentiles, and the line within each box 
indicates the median. Whiskers extend to the highest and lowest values that are not considered outliers (three standard deviations above or below the 
median). Outliers are plotted as individual points. (C) Stress-induced Sko1 binding requires Hog1, but not Hog1 phosphorylation, of Sko1. Sko1 binding 
was measured by ChIP against the HA epitope in wild-type, Ahogl, and slco7- S108A, T1 13A, S126A (which cannot be phosphorylated by Hog1) strains, all 
tagged with three copies of the HA tag at the N' of the endogenous copy of Sko1 . ChIP experiments were performed on cells in log-phase growth in rich 
medium, in the presence or absence of osmotic shock (induced by 5 min in 0.4 M KCI). The data points on the bar graph are the averages of at least three 
biological replicates, and error bars represent the standard deviation of the measurements. 
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recognize specific a'5-regulatory elements at the ORFs that it occupies 
during stress. Hogl interacts directly with RNA Pol II during stress 
(Pascual-Ahuir et al 2006), and is capable of inducing transcription 
when tethered to a promoter (Alepuz et al 2003). Based on these 
observations, we suspect that Hogl recruits RNA Pol II to the set of 
genes with nearby Skol and Hotl binding sites and that Hogl is 
loaded onto these ORFs along with RNA Pol II during transcription 
initiation. 

At the set of genes that show stress-induced recruitment of Skol, 
Hotl, and Hogl to their promoters and/or downstream regions, Hogl 
enrichment in ORFs is largely proportional to RNA Pol II occupancy 
(Figure 5A), consistent with previous reports that Hogl is part of the 
elongating RNA Pol II holoenzyme (Proft et al 2006). This relation- 
ship between gene expression and Hogl presence during elongation is 
specific to stress and does not occur in control conditions (Figure 5B). 
Taken together, these observations suggest that genes marked by the 
presence of Hogl, Skol, and Hotl in their regulatory regions are 
selectively transcribed by a Hogl -RNA Pol II complex during stress. 

DISCUSSION 

Hogl and RNA Pol II were known to colocalize at some stress- 
responsive genes during stress (Pokholok et al 2006; Proft et al 
2006), and Hogl has been shown to associate with components of 
the elongating RNA Pol II holoenzyme under stress conditions (Proft 



et al 2006), but the role of this Hogl -RNA Pol II complex in the 
stress response was not clear. Our results suggest that this complex 
marks and specifies the most highly induced genes during stress. In 
response to stress, cells experience a rapid drop in their overall tran- 
scription rate. This stress-induced global reduction of transcription 
may have two causes: (1) a general effect of the physical changes that 
accompany stress (Proft and Struhl 2004); and (2) a Hogl -dependent 
signal that increases this effect. Given that Hogl interacts with the 
polymerase upon stress (Alepuz et al 2003; Pokholok et al 2006; Proft 
et al 2006), we suggest that the entry of Hogl into the nucleus upon 
stress may serve as this signal. In this model, RNA Pol II either in 
complex with Hogl or modified by Hogl via phosphorylation 
has reduced affinity for non- stress-responsive genes, allowing Hogl 
to briefly hijack some fraction of RNA Pol II to carry out a Hogl- 
specified transcription program. If high levels of nuclear Hogl are 
capable of globally altering transcription, we would expect spurious 
activation of Hogl to be detrimental. In fact, overexpression of Hogl, 
which results in increased levels of nuclear Hogl, is lethal (Wurgler- 
Murphy et al 1997), and genetic perturbations of upstream signaling 
that result in constitutive activation of the Hogl pathway are not 
viable (Maeda et al 1994, 1995). 

The entry of Hogl into the nucleus upon stress allows two Hogl 
cognate transcription factors, Skol and Hotl, to bind with Hogl at 
cis- elements contained in the regulatory regions of stress responsive 
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Figure 5 A stress-specific binding mode of Sko1 targets Hog1 to 
ORFs. (A) In stress, genes bound by a Sko1/Hot1/Hog1 complex (see 
Table S3) show highest levels of RNA Pol II and Hog1 in ORFs. Scatter 
plot shows ChlP-seq of Hog1-HA compared with ChlP-seq of Rpb3 
after 5 min of osmotic shock with 0.4M KCI. The number of reads that 
align to each ORF are normalized to reads per kilobases per million 
(RPKM). Genes that exhibit stress-induced colocalization of Hog1 with 
Sko1 or with both Sko1 and Hot1 in their regulatory region are shown 
in red; all other genes are shown in black. (B) Scatter plot shows ChlP- 
seq of Hog1 -HA compared with ChlP-seq of Rpb3 during rapid growth 
in rich medium. 

genes, targeting the Hog 1 -RNA Pol II complex to these genes. The 
RNA Pol II-Hogl complex transcribes highly expressed, stress- 
induced genes, while avoiding housekeeping genes that were previ- 
ously highly expressed. The selective targeting of this complex suggests 
a mechanism by which a set of genes could be exempted from the 
transient global transcriptional repression that accompanies stress. 

The function of polymerase eviction from highly expressed meta- 
bolic genes, e.g., PDC1, in the early stages of the osmotic stress response 
may be to reallocate transcriptional resources rather than to modify the 
expression level of the RNA Pol II-depleted genes. Depletion of RNA 
Pol II from highly expressed ORFs often does not lead to reduced 
mRNA levels, as many of the transcripts produced by Pol II depleted 
genes are in fact stabilized during the stress response (Miller et al. 2011; 
Romero- Santacreu et al. 2009). Global depletion of RNA Pol II, cou- 
pled with a gene-specific override mechanism, may constitute a strategy 
that allows the cell to rapidly switch its transcriptional program. 
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